Symmetry constraints on phonon dispersion in graphene 
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Taking into account tlie constraints imposed by the lattice symmetry, the phonon dispersion is 
calculated for graphene with interactions between the first and second nearest neighbors in the 
framework of the Born- von Karman model. Analytical expressions are obtained for the out-of-plane 
(bending) modes determined only by two force constants as well as for the in-plane modes with four 
force constants. Values of the force constants are found in fitting to elastic constants and Raman 
frequencies observed in graphite. 
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I. INTRODUCTION 

Since the discovery of graphene (a single atomic 
layer of graphite)^^, main attention has been devoted 
to its electronic properties. More recently, Raman 
spectroscopy^ extends to investigations of graphene. For 
interpretations of the Raman scattering as well as of the 
transport phenomena, the detailed knowledge of the lat- 
tice dynamics and the electron-phonon interactions is 
needed^. 

Several mo dels^'^'^'^'^ ' ^ ° i ^ ^ i ^ ^ have been proposed to 
calculate the phonon dispersion in bulk graphite. Most 
improved ones^iii^ involve many (up to twenty) pa- 
rameters. Recently, detailed measurements and first- 
principles calculations of optical phonon frequencies were 
made for graphit O '^^'^^. They show the qualitative dis- 
agreement with the models^ii^, which include the central 
and angular atomic forces between the first and second 
neighbors in the graphite lattice. 

The passage in the lattice dynamics from graphite to 
graphene and then to nanotubes was examined in the ab 
initio calculationi^ii^, in RefJi using the model^, and in 
Rcf.^'* up to the fourth neighbor with 12 force constants. 
Numerical calculations based on Hiickel's theoryiS and in 
terms of the electron energy were performed in work^S. 
The first-principles calculations^ of dynamical proper- 
ties for graphite and graphene (and also for diamond) 
show that distinctions between the phonon frequencies in 
graphene and relating ones in graphite are negligible in 
comparison with the experimental errors for that frequen- 
cies in graphite. For the highest frequencies, this could 
be intuitively expected because interactions between the 
adjacent layers in graphite are weak. 

Our aim here is to find an analytical description of the 
phonon dispersion in graphene. This can be done within 
the framework of the Born-von-Karman model for the 
honeycomb graphene lattice with interactions only be- 
tween first and second nearest neighbors, but the con- 
straints imposed by the lattice symmetry should be taken 
into account. We show that the out-of-plane (bending) 
and in-plane modes are decoupled from each other. The 
out-of-plane modes are described by two force constants 
determined in fitting to the Raman frequency and small- 



est elastic constant C44. In the narrow wave- vector in- 
terval near the F point, the acoustic out-of-plane mode 
has a linear dispersion with the non-zero sound veloc- 
ity in contrast to the so-called membrane mode which 
demonstrates the quadratic dispersion in the continuous- 
medium approximation^^. We do not pay close attention 
to the agreement of lower frequencies with experiments 
on graphite because their values in graphene are convinc- 
ingly less that ones in graphite. The four force constant 
describing the in-plane modes are found in fitting to their 
Raman frequencies and elastic constant Cn and C12 of 
graphite. The extent of agreement of the present theory 
with experiments corresponds to the comparison level be- 
tween the first-principal calculations^^ and their data (see 
Table 3). 



II. PHONON DYNAMICS IN NEAREST 
NEIGHBOR APPROXIMATION 

The equations of moution in the harmonic approxima- 
tion are written in the well-known form 



E 



(a„ 



Ouf (a„)-c^\"(a„) = 0, (1) 



where the vectors a„ numerate the lattice cells, the su- 
perscripts K, k' note two sublattices A and B, and the 
subscripts i, j = x,y, z take three values corresponding to 
the space coordinates. Since the potential energy is the 
quadratic function of the atomic displacements M^(a„) 
and u^(a„), the force-constant matrix can be taken in 
the symmetric form, $^-^(a„) = <i>^"*(— a„), and its 
Fourier transform, i.e. the dynamical matrix, is a Her- 
mitian matrix. 

Each atom, for instance, Aq (see Fig. 1) has three 
first neighbors in the other sublattice, i. e., B with the 
relative vectors 



Bi = a(l,0), B2,3 



-l,±\/3)/2, 



and six second neighbors in the same sublattice A with 
the relative vectors 



Ai^4 = ±a(0,\/3), 

A3,6 



A2,5 = ±a(-3, V3)/2, 
= Ta(3,%/3)/2, 
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where a = 1.42 A is the carbon-carbon distance. 

For the nearest neighbors (in the B sublattice), the 
dynamical matrix has the form 

'/'^^(q)-ELi'i>4''(BK)exp(zqBJ, (2) 

and for the next neighbors (in the A sublattice) 

^l^^(q) = <i>^.^(Ao) + ELi *^-^(Ak) exp(*qAJ, (3) 

where Aq indices the atom chosen at the center of the 
coordinate system in the A sublattice and the wave vector 
q is taken in units of 1/a. 



FIG. 1: First and second neighbors in the graphene lattice. 

The point group Dg^ of the honeycomb lattice is gen- 
erated by {Cg, di,, CTz}, where <7z is a reflection z —z 
by the plane that contains the graphene layer, Cg is a 
rotation by 7r/3 around the z axis, and is a reflec- 
tion by the xz plane. The transformations of the group 
impose constraints on the dynamical matrix. To obtain 
them, we introduce variables ^,7] = x ±iy transforming 
under the rotation C3 around the z-axis (taken at the Aq 
atom) as follows (^,77) — > (^, 77) exp(±27ri/3). In the ro- 
tation, the atoms change their positions Bi — > B2 — > B3, 
Ai ^ A3 A5, and A2 ^ A4 ^ Ag. Therefore, aU 
the force constants ^'^^(Bf^) with the different k (as well 
as $^^^(Bk)) are equal to one another, but the force con- 
stants with the coincident subscripts ^ or 77 transform as 
covariant variables. For instance. 



AB 



55 (Bi) = $^5^(B2)exp(27r7/3) 
= $ff(B3)exp (-2^7/3). 



(4) 



The relation between $55^ (A^) with the points 
Ai, A3, A5 (and also between A4, A2, Ag) has the same 
form. The constants = $^,^(Bi),72 = ^f^{Ai),a = 



$£'^(Ai) are evidently real. The 



<i>5^,f (Bi), and 7 
constant (3 = $55^ (Bi) is real because the reflection 
{x,y) — > {x,—y) with Bi Bi belongs to the sym- 
metry group. Besides, for the first and second neighbors, 
we have one complex force constant 6 = $55^ (Ai). 

Two force constants $4^(Ao) and ^f^{Ao) for the 
atom Aq can be excluded with the help of conditions 
imposed by invariance with respect to the translations 
of the layer as a whole in the x/z directions. Using the 



equations of motion ^ and Eqs. ([2|), Q, we find this 
stability condition 

(Ao) + 6$4-^(Ai) + 3<i>f,f (Bi) = (5) 

and the similar form for the zz components. 

A. Dispersion of the bending out-of plane modes. 

In the first- and second-neighbor approximation, the 
out-of-plane vibrations uf, uf in the z direction are not 
coupled with the in-plane modes because the force con- 
stants of type ^xz or equals zero due to the reflection 
2: — > —z. The corresponding dynamical matrix for the 
out-of-plane modes has the form 



lAB 



(q) 

0^,^(q)* </)^/(q) 



(6) 



where 



27^[cos {V^Qy) -)- 2cos {3qx/2) cos {\/3qy/2) - 3] - Sa^, 

'/'ff (q) = (7) 
az[exp {iqx) + 2 exp {-iqx/2) cos {\/3qy/2)]. 

The phonon dispersion for the out-of-plane modes is 
found 



o;zo,ZA(q) = V04^(q)±|04^(q) 



(8) 



The equations allow us to express the phonon frequencies 
of the out-of plane branches at the critical points F, K, 
and M in terms of the force constants: 



^zo(r) = \/-6az, tozo, za{K) = ^J-Za^ - 97^, (9) 



t^zo, za(M) = ^(-3Tl)a2-872. 

Expanding Eq. ([8]) in powers of the wave vector q, 
we find the velocity of the acoustic out-of-plane mode 
propagating in the layer 

5, =a[-0.75a^-4.57^]'/' = TWp, (10) 

where we use the well-known formula for the velocity 
of the acoustic z-mode propagating in the x-direction in 
terms of the elastic constant C44 and density p of a hexag- 
onal crystal. Because the interaction between the layers 
in graphite is weak, we can correspond the values of C44 
and p to graphite. Using the values of az and 72 obtained 
in fitting to the exprimental data (see Table 1), we find 
the value of the sound velocity for the out-plane mode 
Sz = 1.57 X 10^ cm/s. 



B. Dispersion of the in-plane modes. 

The dynamical matrix for the in-plane vibrations has 
the form similar to that for the in-plane mode ©, but 
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instead of the functions (j)f^{q} and (f>ff{q} we have to 
substitute respectively the 2x2 matrices 



lAA 



(q) 



(q) 



</>|lf(q) 0f/(q) 



0f,^(q) 0^^(q) 

<,^(q) 0f,f(q) 



(11) 



(12) 



The matrix elements 0^,j^(q) and 4>f^{<l) are obtained 
from 0^/(q) and (f>ff{q), Eqs. ([7]), respectively, with 
substitutions 7 and a instead of 7^ and a^. The off- 
diagonal elements are given by 

^tti^) = (13) 
6[exp{iy/3qy) + 2 cos {3qx/2 + 2tt/3) exp{—iy/3qy/2)] + 
(5*[exp(-i\/3qj^) + 2 cos (3g:r/2 - 27r/3) exp (iV3(7y/2)], 

-/-ff (q) = 

/3[exp (iga;) + 2 exp {-iqx/2) cos (\/3(7.y/2 - 27r/3)]. 

The matrix elements for the B sublattice can be ob- 
tained from that ones for the A sublattice by C2 rotation 
(x, y) — > ~{x, y) of the graphene symmetry group. 

The optical phonon frequencies for the in-plane 
branches at F and K are found 



u> 



in— pi 



1.2 



(r) 



in — pi 
■^1,2 



V— 60:, doublet, 
{K) = v/-3a-97, doublet, 

in— pi 
3.4 



(14) 



{K) ^ yj-ia - 97 ± 3/3. 



Using Eqs. (fTT|) - (fT3| . we find the in-plane mode disper- 
sion in the explicit form for the G ~ K direction. An 
algebraic equation of the forth order have to be solved 
for the M point as well as for points of the general posi- 
tion. 

The in-plane vibrations make a contribution into the 
elastic constants Cn and C12. The corresponding rela- 
tion between the dynamic matrix elements and the elas- 
tic constants can be deduced taking the long-wavelength 
limit (q 0) in the matrices (fTTj) and (fT^ . In this limit, 
separating the acoustic vibrations u'*'^ from the optical 
modes, we obtain the equation of motion in the matrix 
form 



AB 



iBB 



F0^^)/2 

,21 „ac ^ 



U" 



(15) 



where the subscripts and 1 mean that the terms of the 
zero and first order in q should, correspondingly, be kept 
in the matrices pT|) and (|12p but the expansion to the 
second order is used in other terms. We find the matrix 
factor of u'"'' in Eqs. p5)) : 



siq 



- LO 



S2C 



S2q^ 

siq^ - 



(16) 



where 



Sl 



-^7- J(«-/3V«), S2 = ^5-^/3. (17) 



With the help of Eq. (|T6)) . we obtain the velocities of 
longitudinal and transverse acoustic in-plane modes 



SLA = a^/si -I- S2 = \JCi\Ip, 

STA = a^JS1 - S2 = a/(Cii - Ci2)/2/9, 



(18) 



corresponding them to elastic constants and density of 
graphite. The values of force constants from Table 1 
give the sound velocities s la = 2.03 x lO^cm/s, sta = 
1.31 X lO^cm/s. 



TABLE I: Force constants in 10^ cm 



a 


/3 


7 


S 






-3.980 


-1.132 


-0.297 


1.123 


-1.270 


0.204 



TABLE II: Elastic constants (in 10 GPa) calculated (theo) 
and observed-'' 



(exp). 





Cii 


C12 


C44 


theo 


93 


16 


0.56 


exp 


106 ±2 


18 ±2 


0.45 ± .05 



The calculated phonon dispersion is shown in Fig. 2. 
Notice, first, that the sound velocities (for the long waves, 
q — > F) are isotropic in the xy plane as it should be ap- 
propriate for the symmetry of graphene. Second, the in- 
plane LO/TO modes at F, the in-plane LO/LA modes 
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FIG. 2: Calculated phonon dispersion for graphene; the force 
constants, elastic constants, and phonon frequencies in critical 
points are listed in Tables 1, 2, and 3 correspondingly. 
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TABLE III: Phonon frequencies at critical points in cm ^; z and || stand for the out-of-plane and in- plane branches, respectively. 





r [0 0] 


M [l^]7r/3a 


K [01]47r/3V3a 






'^1 '^2 '^3 '^4 '^1 


1 2,3 4 '-*-^l,2 


exp 


1590" 861" 
1583'' 868"= 
1565*' 874" 


1389" 630'* 670" 471" 
1390'' 1323'' 1290'' 625" 625" 451"* 

480" 


1313'' 1184'' 482^* 
1265'' 1194'' 517^* 
1285" 1011" 531" 


theo*" 


1581 


1425 1350 1315 


1300 1220 950 


theo 


1545 873 


1432 1307 1268 598 587 301 


1342 1209 1059 444 



" Reference^^, '' Reference^^, " Reference^, Reference^^, " Reference^^ 



at if, and the out-of-plane ZA/ZO modes at K are dou- 
bly degenerate, because graphene is the non-polar crys- 
tal and the C^v symmetry of these points in the Brillouin 
zone admits the two-fold representation (observation of 
splitting of that modes in graphene would display the 
symmetry braking of the crystal). 

III. FITTING TO EXPERIMENTAL DATA FOR 
GRAPHITE 

Because of the lack of information on graphene, 
we compare the present theory with experiments on 
graphite. In fitting to the experimental data, we must 
keep in mind that the frequencies in graphene for the 
out-of-plane branches should be less than their values in 
graphite, since the atoms are more free to move in the z 
direction in graphene comparatively with graphite. It is 
evident that the adjacent layers in graphite affect more 
intensively the low frequencies. The interaction of the 
adjacent layers can be estimated from the observed split- 
ting at around 130 cm"'^ [see, for instance, Ref.^^] of the 
ZA and ZO' modes which slick together if the interaction 
is absent. 

Thus, we have only two force constants and 7^ to fit 
four Raman frequencies of the out-of-plane modes. We 
see also from Eq. pUj) that these force constants deter- 
mine the velocity Sz of the acoustic out-of-plane mode 
along with the elastic constant C44 . This our result con- 
tradicts to the well-known statement^ that the acoustic 
out-of-plane mode has a quadratic dispersion. That can 
be found only if the condition 7^ = — a^/G of the veloc- 
ity vanishing is satisfied. Therefore, the force constant 
az given in Table 1 is obtained using the experimental 
value of the out-of-plane Raman frequency. The small 
force constant 7^ is found from the value of C44, Table 2, 
taking the experimental error into account. As one can 
see in Table 3, the lower frequencies of the out-of-plane 
mode at other critical points in graphene turn out to be 
less than their values in graphite. The interlayer inter- 
action should diminish that differences. On the other 
hand, the fact that the sound velocity is very sensi- 
tive to the small variation of 72 and becomes complex for 
72 > 0.21 X 10^ cm~^ indicates that graphene is nearly 
unstable with respect to transformation into a phase of 



the lower symmetry group at F. 

For the in-plane modes, we have to fit eight Raman 
frequencies and two elastic constants using four force 
constants. The equations (|17p . and can be 

used as a starting point. Fitting of the in- plane branches 
is insensitive to the imaginary part of the constant 8. 
Therefore, it is taken as a real parameter. Results of fit- 
ting are presented in Fig. 2 and in Tables. Notice, that 
the extent of agreement of the present theory with the 
data obtained for graphite corresponds to the comparison 
level between the first-principal calculations for graphite 
in Refji^ and their experimental data (see Table 3). We 
see only qualitative discrepancy in the subsequence of the 
levels at M: in Fig. 2 the highest level is the LO mode, 
whereas in ReL^S the crossover of the TO and LO modes 
is found out on the F-M line (similar to F-K line), yield- 
ing the TO mode higher at M. We examined versions 
with the crossover. The agreement with experiments is 
not so good in these cases as for one that shown in Fig. 
2 and Tables, but the discrepancies of the order of 50 
cm~^ between the different experiments as well as the 
distinctions between graphene and graphite do not allow 
us to choose the version conclusively. The experiment on 
graphene would clarify this point. 



IV. SUMMARY 



We calculate the phonon dispersion in graphene us- 
ing the Born-von-Karman model with only the first- and 
second-neighbor interactions imposed by the symmetry 
constraints. The bending (out-of-plane) modes are not 
coupled with the in-plane branches and indicate the la- 
tent instability of graphene with respect to transforma- 
tion into a lower-symmetry phase. The Raman frequen- 
cies of these modes are less than the corresponding val- 
ues in graphite. The fitting of the higher in-plane modes 
shows the good agreement of the calculated optical fre- 
quencies as well as elastic constants with experiments. 
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